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Abstract 



We consider simple stochastic climate models, described by slowly time-dependent 
Langevin equations. We show that when the noise intensity is not too large, these 
systems can spend substantial amounts of time in metastable equilibrium, instead of 
adiabatically following the stationary distribution of the frozen system. This behaviour 
can be characterized by describing the location of typical paths, and bounding the 
probability of atypical paths. We illustrate this approach by giving a quantitative 
description of phenomena associated with bistability, for three famous examples of 
simple climate models: Stochastic resonance in an energy balance model describing 
Ice Ages; hysteresis in a box model for the Atlantic thermohaline circulation; and 
bifurcation delay in the case of the Lorenz model for Rayleigh-Benard convection. 
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1 Introduction 

One of the main difficulties of realistic climate models is that they involve a huge number 
of interacting degrees of freedom, on a wide range of time and length scales. In order 
to be able to control these models analytically, or at least numerically, it is necessary 
to simplify them by eliminating the less relevant degrees of freedom (e. g. high-frequency 
or short-wavelength modes). A possible way to do this is to average the equations of 
motion over all fast degrees of freedom, a rather drastic approximation. As proposed by 
Hasselmann |]20[| (see also (2|) , a more realistic approximation is obtained by modeling the 
effect of fast degrees of freedom by noise. 

In a number of cases, it is appropriate to distinguish between three rather than two time 
scales: Fast degrees of freedom (e.g. the "weather"), which are modeled by a stochastic 
process; intermediate "dominant modes" (e. g. the average temperature of the atmosphere) 
whose dynamics we want to predict; and slow degrees of freedom (e. g. the mean insolation 
depending on the eccentricity of the Earth's orbit), which evolve on very long time scales 
of several centuries or millennia, and can be viewed as an external forcing. Such a system 
can often be modeled by a slowly time-dependent Langevin equation 



dx t = f(x t , et) dt + aG{et) dW t 




where the adiabatic parameter e and the noise intensity a are small parameters, Wj is a 
standard vector-valued Wiener process (describing white noise) and G is a matrix. 



Our aim in this paper is to describe the effect of the noise term on the dynamics of ( |l.l|) , 
assuming the dynamics without noise is known. For this purpose, we will concentrate on 
bistable systems, which frequently occur in simple climate models: For instance, in models 
for the major Ice Ages, where the two possible stable equilibria correspond to warm and 
cold climate 0, or in models of the Atlantic thermohaline circulation |34|, |3(J. Noise 
may enable transitions between the two stable states, which would be impossible in the 
deterministic case, and our main concern will be to quantify this effect. 

The method used to study the stochastic differential equation (SDE) ([O]) will depend 
on the time scale we are interested in. Let us first illustrate this on a static one-dimensional 
example, namely the overdamped motion in a symmetric double- well potential: 

dx t = —^-V(x t )dt + adW t , V(x) = \bx i -\ax 2 . (1.2) 
ox 4 2 

where a and b are positive constants. The potential has two wells at ±\/a/b, separated by 
a barrier of height H = a 2 /(4b). A first possibility to analyse this equation is to compute 
the probability density p(x,t) of xt- It obeys the Fokker-Planck equation 



d . , d 



dV 

— (x)p(x,t) 



r 2 Q2 



■ Y7J^'' UJ) - (L;!! 



which admits in particular the stationary solution 



Po(x) = ±e- 2V ^\ (1.4) 

where N is the normalization. At equilibrium, there is equal probability to find xt in 
either potential well, and for weak noise it is unlikely to observe xt anywhere else than in 
a neighbourhood of order a of one of the wells. 

Assume now that the initial distribution xq is concentrated at the bottom yf a/b of the 
right-hand potential well. Then it may take quite a long time for the density to approach its 



asymptotic value (1.4). A possible way to investigate this problem relies on spectral theory. 



Denote the right-hand side of (F3) as £p(x,t), where £ is a linear differential operator. 
The stationary density ( |1.4| ) is an eigenfunction of L with eigenvalue 0. We may assume 
that C has eigenvalues • • • < A& < • • • < A2 < Ai < 0, c. f. pi], Section 6.7]. Decomposing 
p(x,t) on a basis of eigenfunctions of C, we see that p approaches the stationary solution 
in a characteristic time of order l/|Ai|. 

There exists, however, a much more precise description of the process xt than by its 
probability density. Recall that for almost every realization Wt(oj) of the Brownian motion, 
the sample path 1 1— ► xt(u>) is continuous. Instead of computing the time needed for p(x, t) 
to relax to po(x), we can consider the random variable 

t(lo) = inf{t > 0: x t (u) < 0}, (1.5) 

describing the first time at which the path xt crosses the saddle (one could as well consider 
the first time the bottom of the left-hand well is reached) . The distribution of r is asymp- 
totically exponential, with expectation behaving in the weak-noise limit like Kramers' 
time 

^Kramers — 6 ^ . (1 .6) 

A mathematical theory allowing to estimate first-exit times for general n-dimensional 
systems (with a drift term not necessarily deriving from a potential) has been developed 
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by Freidlin and Wentzell [jlgl . In specific situations, more precise results are available, for 
instance subexponential corrections to the asymptotic expression (|1.6|), see (3], Even 
the limiting behaviour of the distribution of the first-exit time from a neighbourhood 
of a unique stable equilibrium point has been obtained The first-exit time from a 

neighbourhood of a saddle has been considered by Kifer in the seminal paper p5|]. 

If the noise intensity a is small (compared to the square root of the barrier height), 
then the time needed to overcome the potential barrier is extremely long, and the time 
required to relax to the stationary distribution po(x) is even longer. In fact, on time 
scales shorter than Kramers' time, solutions of ( |1.2Q starting in one potential well will 
hardly feel the second potential well. As we will see in Section g xt is well approximated 
by an Ornstein-Uhlenbeck process, describing the overdamped motion of a particle in a 
potential of constant curvature c = 2a. The Ornstein-Uhlenbeck process relaxes to a 
stationary Gaussian process with variance a 2 /(2c) in a characteristic time 

Trelax = ~- (1-7) 

c 

Thus for < t <. ^reiaxi the behaviour of xt is transient; for r re iax ramors) Xt 

is close to a stationary Ornstein-Uhlenbeck process with variance a 2 /(2c); and only for 
ramers will the distribution of x% approach the bimodal stationary solution (|1.4|) . 



This phenomenon, where a process seems stationary for a long time before ultimately 
relaxing to a new (possibly stationary) state, is known as metastability. It is all the more 
remarkable in an asymmetric double-well potential: then a process starting at the bottom 
of the shallow well will first relax to a metastable distribution concentrated in the shallow 
well, which is radically different from the stationary distribution having most of its mass 
concentrated in the deeper well. 

A different approach, based on the concept of random attractors (see ]12|. [3l| , [IJ), gives 
complementary information on the long-time regime. In particular, in |l3| it is proved 
that for arbitrarily weak noise, paths of (|l]^) with different initial conditions but same 
realization of noise almost surely converge to a random point. The time needed for this 
convergence, however, diverges rapidly in the limit a —* 0, because paths starting in dif- 
ferent potential wells are unlikely to overcome the potential barrier and start approaching 
each other before Kramers' time. 

We now turn to situations in which the potential varies slowly in time. For simplicity, 
we will consider the family of Ginzburg-Landau potentials 

1 1 

U(x,A,/x) = -x 4 — ^x 2 — Ax, (1.8) 

and let either A or /i vary in time, with low speed e. For instance, A or ^ may depend 
periodically on time, with low frequency 2-ke. The potential V has two wells if 27A 2 < 4// 3 
and one well if 27A 2 > 4^/ 3 , and when A or fi are varied, the number of wells may change. 
Crossing one of the curves 27A 2 = 4/u 3 , \x > 0, corresponds to a saddle-node bifurcation, 
and crossing the point A = fi = corresponds to a pitchfork bifurcation. 

The slow time-dependence introduces a new time scale Tf orcing = 1/e. Since curvature 



and barrier height are no longer constant, we replace the definitions (|1.6| ) and (1.7) by 

T (max) _ 2H max /<7 2 „ J T (min) _ 1 /, q n 

^Kramers — e dnu "'relax — „ ' V 1 -^ 

trnax 

where H max denotes the maximal barrier height during one period, and c max denotes the 
maximal curvature at the bottom of a potential well. Here we are interested in the regime 

T->( min ) s~ rp „ T (max) 1f) s 

± relax ^ 1 forcing ^ -L Kramers' K 1 - 1 ^) 
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which means that the process has time to reach a metastable "equilibrium" during one 
period, but not the bimodal stationary distribution. Mathematically, we thus assume that 
£ <C c max and a 2 <C 2iJ max /|loge|. We allow, however, the minimal curvature and barrier 
height to become small, or even to vanish. 

For time-dependent potentials, the Fokker-Planck equation ( |1.3| ) is even harder to solve 
(and in fact, it does not admit a stationary solution). Moreover, random attractors are not 
straightforward to define in this time-dependent setting. We believe that the dynamics 
on time scales shorter than T^^ ers is discussed best via an understanding of "typical" 
paths. The idea is to show that the vast majority of paths remain concentrated in small 
space-time sets, whose shape and size depend on the potential and the noise intensity. 
These sets are typically located in a neighbourhood of the potential wells, but under some 
conditions paths may also switch potential wells. There are thus two problems to solve: 
first characterize the sets in which typical paths live, and then estimate the probability of 
atypical paths. It turns out that these properties have universal characteristics, depending 
only on qualitative properties of the potential, especially its bifurcation points. 

We start, in Section by discussing the simplest situation, which occurs when the 
initial condition of the process lies in the basin of attraction of a stable equilibrium branch. 
For sufficiently small noise intensity, the majority of paths remain concentrated for a long 
time in a neighbourhood of the equilibrium branch. We determine the shape of this 
neighbourhood and outline how coloured noise can decrease the spreading of paths. 

Section || is devoted to the phenomenon of stochastic resonance. We first recall the 
energy-budget model introduced in j| to give a possible explanation for the close-to- 
periodic appearance of the major Ice Ages. This model is equivalent to the overdamped 
motion of a particle in a modulated double-well potential, where the driving amplitude 
is too small to allow for transitions between wells in the absence of noise. Turning to 
the description of typical paths, we find a threshold value for the noise intensity below 
which the paths remain in one well, while above threshold, they switch back and forth 
between wells twice per period. The switching events occur close to the instants of minimal 
barrier height. Several important quantities have a power-law dependence on the small 
parameters, in particular the critical noise intensity, the width of transition windows, and 
the exponent controlling the exponential decay of the probability of atypical paths. 

In Section f|, we start by discussing a variant 11] of Stommel's box model [gj] of 
the Atlantic thermohaline circulation. Assuming slow changes in the typical weather, 
this model also reduces to the motion in a modulated double-well potential, where the 
modulation depends on the freshwater flux. If the amplitude of the modulation exceeds a 
threshold, the potential barrier vanishes twice per period, so that the deterministic motion 
displays hysteresis. Additive noise influences the shape of hysteresis cycles, and may even 
create macroscopic cycles for subthreshold modulation amplitude. We characterize the 
distribution of the random freshwater flux causing the system to switch from one stable 
state to the other one. 

Finally, in Section [|, we consider the Lorenz model for Rayleigh-Benard convection 
with slowly increasing heating. In the deterministic case, convection rolls appear only some 
time after the steady state looses stability in a pitchfork bifurcation. This bifurcation delay 
is significantly decreased by additive noise, as soon as its intensity is not exponentially 
small. 
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2 Near stable equilibria 



Let us start by investigating Equation ( |1.1| ) in the one-dimensional case, i.e., when xt, 
Wt and G(t) = g{t) are scalar. Since we are interested in the dynamics on the time scale 
^forcing = V e > we rescale time by a factor e, which results in the SDE 

dx t = -f(xt,t)dt + ^ F g(t)dW t . (2.1) 

e ye 

The factor l/y/e is due to the diffusive nature of the Brownian motion. 

In this section, we will consider the dynamics near a stable equilibrium branch of /, 
i. e., a curve x*(t) such that 

df 

f(x*(t),t) = and a*(t) = ^(z*(t),t) < -a (2.2) 

for all t, where ao is a positive constant. In the one-dimensional case, / always derives 
from a potential V, and —a*(t) represents the curvature at the bottom x*{t) of a potential 
well. 



In the deterministic case a = 0, solutions of (2.1) track the equilibrium branch x*(t) 



adiabatically. In fact, Tihonov's Theorem [36, |19 
particular solution xf et with an asymptotic expansion of the form 



asserts that for a = 0, (G.lf) admits a 



x*(t) + e^ + 0(e z ). (2.3) 

Since a*(t) is negative, xf et lies a little bit to the left of x*(t) if x*(t) moves to the right, 
and vice versa. The adiabatic solution xf et attracts nearby solutions exponentially fast in 
t/e. 

Consider now the SDE (^^) with positive noise intensity. For the sake of brevity, we 
assume that g is positive and bounded away from zero. In a nutshell, our main result can 
be formulated as follows: Up to Kramer's time, paths starting near Xg et are concentrated 
in a neighbourhood of order ag(t)/ \/\a*(t)\ of the deterministic solution with the same 
initial condition, as shown in Figure |l[ Larger noise intensities and smaller curvatures 
thus lead to a larger spreading of paths. This result holds as long as the spreading is 
smaller than the distance between x*(t) and the nearest unstable equilibrium (i.e., the 
nearest saddle of the potential). 

To make this claim mathematically precise, we need a few definitions. For simplicity, 
we discuss first the particular case xq = XQ 6t . We use the notations 

df f* 
a (t) = -^-(xf e \t), a(t,s) = J a(u)du and a(t) = a(t, 0). (2.4) 



Note that by ( |2.3| ), a(t) = a*(i) + 0{e) is negative for sufficiently small e. The main 
idea is that xt — xf et is well approximated by a generalized Ornstein-Uhlenbeck process, 
with time-dependent damping a(t)/e and diffusion coefficient og{t)/yfe. This process is 
obtained by linearizing the SDE ( |2.1| ) around xf et , and has variance 



T 2 rt 



v(t) = — e 2a{t ' s)/e g(s) 2 ds. (2.5) 
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Figure 1. A sample path of the SDE (2.1), for f(x,t) — a*(t)(x — x*(t)) deriving from 
a quadratic single- well potential, and g(t) = 1. The potential well is located at x*(t) = 
sin(27rf), and has curvature —a*(t) — 4 — 2sin(47rf). Parameter values are e = 0.04 and 
a = 0.025. After a short transient motion, the deterministic solution xf ct tracks x*(t) at 
a distance of order e. The path xt is likely to stay in the shaded set B(h) (shown here for 
h = 3), which is centred at a;^ et and has time-dependent width of order ha/ y/\a*(t)\. 



The function v(t) solves the ordinary differential equation (ODE) ev = 2a(t)v + a 2 g(t) 2 . 
In analogy with ( |2.3|) , this equation also admits a particular solution v(t) satisfying 

m = ^y l [g(tf + 0(e)], (2.6) 

and since a(t) ^ — a$t for t ^ 0, the variance v(t) = v(t) — v{fS) e 2a ^/ £ approaches v(t) 
exponentially fast. We now introduce the set 

B{h) = {(x,t): \x-x^\ < hyftXfi}, (2.7) 

which depends on a real parameter h > 0. The strip B{h) is centred in the adiabatic 
solution xf et tracking the bottom of the potential well, and has time-dependent width 
hag{t)/^2\a*(t)\[l + 0(e)]. To lowest order in e and a, B(h) coincides with the points in 
the potential well for which V(x,t) — V(xf et ,t) is smaller than (^hag(t)) 2 . 

The main result is that for h 3> 1, paths {x s } s ^,q are unlikely to leave the set B(h) 
before Kramers' time. Equivalently, the first-exit time 

T B{h) =mf{t>0: (x t ,t) <£ B(h)} (2.8) 

is unlikely to be smaller than Xkramers- Indeed, one can prove the following estimate (see 
H Theorem 2.4] and [Jsj, Theorem 2.2]). There is a constant ho > such that 

F{r B{h) >t}^C(t,e)e- Kh2 (2.9) 
holds for all t > and all h ^ ho /a, where 

C(t,e) = ^r^- + 2 and k = - - 0(e) - O(ah). (2.10) 

The exponential term e~ Kh2 in ( |2.9| ) is independent of time, and becomes small as soon 
as h S> 1. The constant ho depends on / and is the smaller the smaller ao is: The flatter 
the well, the more restrictive the condition h ^ ho jo becomes. The prefactor C(t,e), 
which grows as time increases (and is certainly not optimal) only leads to subexponential 
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corrections on the time scale Tf orc i ng . Some time dependence of the prefactor is to be 
expected, as it reflects the fact that occasionally a path will make an unusually large 
excursion, and the longer we wait the more excursions we will observe. The prefactor 
also depends on e. A factor 1/e is due to the fact that we are working on the time scale 
^forcing, while the actual factor of 1/e 2 allows us to obtain the best possible exponent. 
Choosing k slightly smaller allows to replace e 2 by the more natural e in the definition 
of C(t,e). Thus we find that paths are unlikely to leave B(h) before time t, provided 
hi/ a 2 > h 2 » log C(t,e). 

The same results hold if xt does not start on the adiabatic solution xf et , but in some 
deterministic xq sufficiently close to it. Then xf et has to be replaced in (|2.4| ) and (|2.7|) 
by the solution xf et of the deterministic equation with initial condition xq. We still have 
that a(t) is negative (and bounded away from zero), but note that fl2.6j ) may not hold for 
very small t, when x^ et has not yet approached xf et . 

If the potential V grows at least quadratically for large \x\, one can deduce from ( |2.9| ) 
that the moments of \xt — xf et \ are bounded by those of a centred Gaussian distribution 
with variance of order v(t), for times small compared to Kramers' time |8], Corollary 2.4], 
even if V has other potential wells than the one at x*(t). Assume for instance that V has 
two potential wells, with the shallower one at x*(t). Then the system is in metastable 
"equilibrium" for an exponentially long time span during which the existence of the deeper 
well is not felt. 

Similar statements are valid in the multidimensional case (in which / does not neces- 
sarily derive from a potential). Let x*(t) be an equilibrium branch of /, and denote by 
A*(t) the Jacobian matrix of / at x*(t). We assume that the eigenvalues of A*(t) have 
real parts smaller than some negative constant —clq for all times, so that x*(t) is asymp- 
totically stable. In the deterministic case a = 0, Tihonov's theorem shows the existence 
of an adiabatic solution 

xf et = x *(t) + eA*®- 1 **® + 0(e 2 ), (2.11) 

which attracts nearby orbits exponentially fast. Let A(t) be the Jacobian matrix of / at 
xf et . It satisfies A(t) = A*(t) + 0{e). The solution of the SDE flD]) linearized at xf ct has 
a Gaussian distribution, with covariance matrix 

a 2 /"* 

X(t) = — U{t,s)G(s)G(s) T U{t,s) T ds, (2.12) 
e Jo 

where U(t, s) is the fundamental solution of ey = A(t)y with initial condition U(s, s) = 1. 
To keep the presentation simple, we will assume that the smallest eigenvalue of G(s)G(s) T 
is bounded away from zero and the largest one is bounded above. Note that X(t) obeys 
the ODE eX = AX + XA T + a 2 GG T , and approaches exponentially fast a matrix X(t) 
which satisfies 

X(t) = X (t) + 0(e), where AXq + XqA t = —a 2 GG T . (2.13) 

Given a deterministic solution xf et , the definition of the set B{h) reads now 

B{h) = {{x,t): {x-xf*) T X(t)-\x-a$*) < h 2 }, (2.14) 

and (|2.9| ) generalizes to the following statement (see || Theorem 6.1] for a discussion; the 
proof will be given in (5|): There is a constant ho > such that for all h ^ ho/<J and all 
«€ (0,1/2), 

Hr B{h) >t}^ C(t,e)e- fth ^ 1 -°^-°^\ (2.15) 
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where 



C(i, e )= - + 1 



1 



n/2 



(2.16) 

n being the dimension of x. Paths are thus concentrated, up to a given time t, in sets of the 
form B(h), which have an ellipsoidal cross-section defined by X(t). Again the parameter 
h must satisfy hi/a 2 > h 2 » logC(i, e). 

This result can be used, in particular, to understand the effect of coloured noise. 
Assume for instance that the one-dimensional system 

dx t = f(x t ,et)dt + g(et)dZ t (2.17) 
is not driven by white noise, but by an Ornstein-Uhlenbeck process Zt obeying the SDE 



dZ t = -'jZtdt + adWt, 



(2.18) 



The equations ( |2.17| ) and ( 2.1S| ) can be rewritten, on the time scale 1/e, as a two- 



dimensional system of the form ( |1.1| ) for (xt,Zt) 
librium branch x*{t) with linearization a*(t) ^ 
asymptotic covariance matrix ( |2.13| ) is given by 



We assume that / has a stable equi- 
-ao < 0. To leading order in e, the 



X (t) 



(7 



g(t) 



\ 



2( 7 + |o*(t)|) 
9(t) 



2( 7 + |o*(t)|) 
1 

27" 



(2.19) 



V2( 7 +|a*(t)|) 

The conditions on GG T mentioned above can be relaxed (c. f. |8], Theorem 6.1]), so 
that ( 2.15 ) is applicable. We find in particular that the path {xt}t^o is concentrated 
in a strip of width proportional to erg if)/ -v/7 + |a*(i)|, centred around xf ct . Hence larger 
"noise colour" 7 yields a smaller spreading of the paths, in the same way as if the curvature 
of the potential were increased by 7. 



3 Stochastic resonance 



In the previous section, we have seen that on a certain time scale, paths typically remain in 
metastable equilibrium. With overwhelming probability, they are concentrated in a strip 
of order crg(t)/ y/\a*(t)\ near the bottom of a potential well with curvature |a*(i)|. This 
roughly holds as long as the strip does not extend to the nearest saddle of the potential. 
New phenomena may occur when this hypothesis is violated, either because the noise 
coefficient ag(t) becomes too large, or because the curvature or the distance to the saddle 
become too small. Then paths may overcome the potential barrier and reach another 
potential well. This mechanism has various interesting consequences, one of them being 
the effect called stochastic resonance. 

Stochastic resonance (SR) was initially introduced as a possible explanation for the 
close-to-periodic appearance of the major Ice Ages [f|. While this explanation remains 
controversial, SR has been detected in several other physical and biological systems, see 
for instance |2^, |39| for a review. 

The original model in Q is based on an energy balance of the Earth in integrated form. 
The evolution of the mean surface temperature T is described by the differential equation 

j rji 

c— = Q (1 + A cos ujt){l - a(T)) - E(T). (3.1) 
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Here the term R^ n = Q(l + Acostot) is the incoming solar radiation, where Q denotes 
the solar constant, and the periodic term models the effect of the Earth's varying orbital 
eccentricity. The amplitude A of this modulation is very small, of the order 5 x 10~ 4 , 
while its period 2tt/uj equals 92 000 years. The outgoing radiation R out = a(T)R- m + E(T) 
depends on the albedo a(T) of the Earth and its emissivity. c denotes the heat capacity. 

To account for the existence of two stable climate states (warm climate and Ice Age), 
the right-hand side of ( |3.1| ) should have two stable and one unstable equilibrium points. 
The authors of [|| postulate that 

7 (T, = - a (T„ - 1 = ,(l - |) (l - |) (l - |), (3.2) 

where T\ = 278.6 K and T 3 = 288.6 K are the representative temperatures of the two 
stable states, and T 2 = 283.3 K represents the unstable state. Since E(T) ~ T 4 varies 
little on this range, the problem can be further simplified by neglecting the T-dependence 
of E(T) ~ (E). Equation (|]|) becomes 



dr _ (E) 



S[ 1 -^Vl-^Vl-£Vl+^coso;t)+Acosa;t 



The parameter (3 is related to the relaxation time t ~ 8 years of the system via 



(3.3) 



(3.4) 



Let us now transform this system to a dimensionless form. We do this in two steps: First 
we scale time by a factor uj/2tt, so that in the new variables, the system has period 1. 
Then we introduce the variable x = (T - T 2 )/AT, where AT = (T 3 - T t )/2 = 5 K. The 
resulting system is 

dx 1 

— = - \—x(x — xi)(x — X3XI + A cos 27rt) + K cos 27rtl , (3-5) 
at e 

where x x = (T x - T 2 )/AT ~ -0.94 and x 3 = (T 3 - T 2 )/AT ~ 1.06. The adiabatic 
parameter e is given by 

LOT 2(T 3 - T 2 ) 4 , , 

e = \^ ^ 1.8 x 10~ 4 . 3.6) 

2vr AT y J 

This confirms that we are in the adiabatic regime. Using the value {E)/c = 8.77 x 
10~ 3 /4000 Ks" 1 from 0], we find a driving amplitude 

_ A JYT2J3 

(3 (AT)3 - U - i2 ' [6J) 

The term in brackets in (|3.5[ ) derives from a double-well potential, which is almost of the 
Ginzburg-Landau type ([Qj|). If we set, for simplicity, x\ = —1 and x 3 = 1, and neglect the 
term ^4cos27ri, then we obtain indeed a force deriving from the potential ( |1.8| ), with fi = 1 
and A = i^cos2-7rt. This potential has two wells if and only if |A| < A c = 2/3\/3 — 0.38, 
and thus the amplitude K of the forcing is too small to enable transitions between the 
potential wells. Note, however, that although A is very small, K is not negligible compared 
to A r . 



9 



u u w 

Figure 2. The potential V(x, t) = jx 4 — ^x 2 — K cos(2-Kt)x, from which derives the drift 
term in (3.8). For cos(27rf) = 0, the potential is symmetric (middle), for integer times, the 
left-hand well approaches the saddle (right), while for half-integer times, the right-hand 
well approaches the saddle (left). If the amplitude K is smaller than the threshold A c , 
there is always a potential barrier, which an overdamped particle cannot overcome in the 
deterministic case. Sufficiently strong noise, however, helps the particle to switch from the 
shallower to the deeper well. This effect is the stronger the lower the barrier is, so that 
switching typically occurs close to the instants of minimal barrier height. 



The main new idea in [BJ] is that if one models the effect of the "weather" by an additive 
noise term, then transitions between potential wells not only become possible but, due to 
the periodic forcing, these transitions will be more likely at some times than at others, so 
that the evolution of T can be close to periodic. We will illustrate this on the model SDE 

dx t = -[x t - x z t + K cos 2vrt] dt + dW t . (3.8) 

e ye 

However, the results i n fj\ apply to a more general class of periodically forced double- well 
potentials, including (|3.5|) . 

Various characterizations of the effect of noise on the dynamics of ([T^), and various 
measures of periodicity have been proposed. A widespread approach uses the signal-to- 
noise ratio, a property of the power spectrum of xt, which shows peaks near multiples of 
the driving frequency [16, p?|, 24]. For small driving amplitudes K, the signal-to-noise ratio 



behaves like e~ H l a /a 4 , where H is the height of the potential barrier in the absence of 
periodic driving (i. e., for K = 0). The signal's "periodicity" is thus optimal for a 1 = H/2. 
A different approach is used in |jl7|| , where the L p -distance between sample paths and a 
periodic limiting function is shown to converge to zero in probability as a — > 0. This result 
requires e to be of order Q- 2H / a ; which implies exponentially long forcing periods. 

We examine here a different regime, in which the forcing amplitude K is not necessarily 
a small parameter, but may approach A c . In this way, transitions become possible for 
values of e which are not exponentially small. The potential barrier is lowest at integer 
and half-integer times. At integer times, the left-hand well approaches the saddle, while 
at half-integer times, the right-hand well approaches the saddle, c. f. Figure 0. 

The minimal values -ffmiiu c m j n and <5 m ; n of the barrier height, the curvature at the 
bottom of the wells, and the distance between the bottom of one of the wells and the 
saddle can be expressed as functions of a parameter ao = A c — K. For small ao, they 

3/2 1/2 1/2 3/2 

behave like -ffmin^oo) x a , c m i n (ao) x a and 5 m in(ao) ^ a o (meaning c_a ^ 
-f^min(ao) ^ c + a o^ f° r some positive constants c± independent of ao, and so on). 

Intuitively, our results from Section || indicate that the maximal spreading of paths 
is of order c/c m i n (ao)^^, provided this value is smaller than <5 m i n (a ), i.e., provided a <C 
a^ 4 . Assume for instance that we start at time 1/4 (when the potential is symmetric) 
near the right-hand potential well. We call transition probability the probability Ptrans 
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of having reached the left-hand potential well by time 3/4, after passing through the 
configuration with the shallowest right-hand well. Extrapolating ( |2.9| ) with h of the order 
KinclliJcJ x H^Ja, we find 

P tra „s £ ^ e"" "* ^ ^ = ^ H ^ for a < af . (3.9) 

Note the similarity with Kramers' time for the potential frozen at the moment of minimal 
barrier height. 

A bound of this form can indeed be proved, but ( |3.9|) turns out to be a little bit 
too pessimistic for very small do- This is a rather subtle dynamical effect, related to the 
behaviour of the deterministic system. Recall that the set 13(h) in ( |2.7| ) is defined via 
the linearization at the adiabatic solution x^ et , not at the bottom x*(t) of the potential 
well. This distinction is irrelevant as long as the minimal curvature remains of order one, 
but not when it is a small parameter. In that case, the asymptotic expansion ( |2.3j ) does 
not necessarily converge. Using methods from singular perturbation theory []n|, one can 
show that xf et never approaches the saddle closer than a distance of order y/e, so that the 
curvature at xf et never becomes smaller than a quantity of order y/e, even if tio < e. As a 
consequence, for oq < e, the system behaves as if there were an effective potential barrier 
of height e 3//2 . 

In fact, one can prove the following bound (see |7], Theorem 2.6] and Theorem 3.1]): 
There exist constants C, k > such that 

Ptrans < J e"^/- 2 for (T ^ <7 C = (do V ef' 4 , (3.10) 

where a V b denotes the maximum of two real numbers a and b. In addition, paths remain 
concentrated in a set B(h) of the form fl2.7| ). Examining the behaviour of the integral fl2,5|) , 
one can show that the width of B(h) behaves, near t = 1/2, like ha/{\t-l/2\ x ' 2 Vcr c 1/3 ). The 
various exponents entering these relations do not depend on the details of the potential, 
but only on some qualitative properties of the "avoided bifurcation" , and can be deduced 
geometrically from a Newton polygon [jTo|1 . 

What happens when a exceeds the threshold value cr c ? Away from half-integer times, 
the right-hand well may still be sufficiently deep to confine the paths. However, there are 
time intervals near half-integer t during which it becomes possible to overcome the barrier. 
Near t = 1/2, the curvature c(t) at xf et and the distance between xf et and the saddle both 
behave like \t — 1/2 1 V at^ . Transitions thus become possible for \t — 1/2] ^ a 2 / 3 . 

During this time interval, the process xt makes a certain number of attempts to over- 
come the barrier. If the saddle is reached, xt has roughly equal probability to fall back into 
the right-hand well, in which case it will make further attempts to cross the barrier, or to 
fall into the deeper left-hand well, where it is likely to stay during the next half-period. 
One can show that the typical time for each excursion is of order e/c(t). Although the 
different attempts are not independent, the probability not to reach the left-hand well 
during the transition window \t — 1/2 1 ^ <r 2 / 3 behaves roughly like (1/2) N , where is the 
maximal number of possible excursions. 

These arguments can be used to show (see 0, Theorem 2.7] and @, Theorem 3.1]) that 
there exist constants C, k > such that 

Ptrans > 1 -Cz~™^l^^ for O > * c . (3.11) 

The factor a 4 / 3 is proportional to the integral of c(t) over the transition window, and the 
factor | logo" | takes into account the time needed to travel from the saddle to the left-hand 
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FIGURE 3. Sample paths of the SDE (|3~8|) for e — a — 0.005, and a = 0.02 (upper 



picture) and a = 0.14 (lower picture). Full curves represent the location of potential wells, 
the broken curve represents the saddle. For weak noise, the path xt is likely to stay in the 
shaded set B(h), centred at the deterministic solution tracking the right-hand well. The 
maximal width of 13(h) is of order ha/(ao V e) 1 / 4 and is reached at half-integer times. For 
strong noise, typical paths stay in the shaded set which switches back and forth between 
the wells at integer and half- integer times. The width of the vertical strips is of order 
er 2 / 3 . The "bumps" are due to the fact that one of the wells becomes very flat during the 
transition window so that paths might also make excursions away from the saddle. 



well. Amplification by SR is thus optimal for noise intensities just above the threshold a c , 
because stronger noise intensities will gradually blur the signal. 

In the large-noise regime a ^ a c , the vast majority of paths stay in a strip switching 
back and forth between potential wells each time the barrier height becomes minimal, as 
shown in Figure ||[ 

Paths spend approximately half the time (for 1/4 < t < 1/2, 3/4 < t < 1, and so on) 
in metastable equilibrium in the shallower potential well. This differs from the quasistatic 
picture, when the driving period is larger than the maximal Kramers time, and paths 
spend most of the time in the deeper potential well with occasional excursions to the 
shallower one. 

While the details of the transition process depend on the potential, the exponents 
in ( |3.10 ) and ( |3.11 ) depend only on qualitative properties of the avoided bifurcation. 
Other exponents arise, for instance, if V is a symmetric potential with modulated barrier 
height of the form (|l.8| ) with A = and pt(t) = do + 1 — cos 2vrt, c.f. g Theorem 3.2]. 
Here an additional feature can be observed: For sufficiently strong noise, the process is 
likely to reach the saddle during a certain transition window, but due to symmetry, it 
has about equal probability to be in either of the wells when transitions become unlikely 
again. Observing the process for several periods, we see that near the instants of minimal 
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barrier height, the process chooses randomly between potential wells, with probability 
exponentially close to 1/2 for choosing either. 

One can also consider the effect of coloured noise on SR. If the system is driven by 
an Ornstein-Uhlenbeck process with damping 7, the typical spreading of paths will be 
smaller, making transitions more difficult. One can show that transitions only become 
likely above a threshold noise intensity <r c , given by 

al = (a Ve)( 7 V(a Ve) 1 / 2 ). (3.12) 

If 7 < (aflVe) 1 / 2 , we recover the white-noise result, but for larger 7, the threshold grows 
linearly with 7, namely like (ao V 5)7. 

It is, of course, not easy to decide whether the observed periodicity in the appearance 
of Ice Ages can be explained by a simple, one-dimensional SDE of the form ( |3.8| ). Our 
results show, however, that in order to match the observations, the noise intensity should 
lie in a relatively narrow interval. Too weak noise will not allow regular transitions between 
stable states, while too strong noise increases the width of the transition windows so much 
that although switching does occur, no periodicity can be observed. 



4 Hysteresis 

The glacial cycle is not the only important bistable system in climate physics. Another 
wellknown example is the Atlantic thermohaline circulation. At present time, the Gulf 
Stream transports enormous amounts of heat from the Tropics as far north as the Barents 
Sea, causing the current mild climate in Western Europe. It is believed, however, that this 
has not always been the case in the past, and that during long time spans, the thermohaline 
circulation was locked in a stable state with far less heat transported to the North (see 
for instance |30tl ). 

A simple model for oceanic circulation showing bistability is Stommel's box model |34|] , 
where the ocean is represented by two boxes, a low-latitude box with temperature T\ and 
salinity Si, and a high-latitude box with temperature T2 and salinity S^. Here we will 



follow the presentation in [11], where the intrinsic dynamics of salinity and of temperature 
are not modeled in the same way. The differences AT = T\ — T2 and AS = S± — S2 are 
assumed to evolve according to the equations 

^-AT = --(AT-6)-Q(Ap)AT (4.1) 
dt r r 

^ AS = l^ls - Q(Ap)AS. (4.2) 

Here r r is the relaxation time of AT to its reference value 9, Sq is a reference salinity, and H 
is the depth of the model ocean. F(t) is the freshwater flux, modeling imbalances between 
evaporation (which dominates at low latitudes) and precipitation (which dominates at 
high latitudes). The dynamics of AT and A5 are coupled via the density difference Ap, 
approximated by the linearized equation of state 

Ap = a s AS - a T AT, (4.3) 

which induces an exchange of mass Q(Ap) between the boxes. We will use here Cessi's 
model for Q, 

Q(Ap) = - + |Ap 2 , (4.4) 
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where t<i is the diffusion time scale, q the Poiseuille transport coefficient and V the volume 
of the box. Stommel uses a different relation, with Ap 2 replaced by |Ap|, but we will not 
make this choice here because it leads to a singularity. 

Using the dimensionless variables y = asAS/(aT0), z = AT/9 and rescaling time by 
a factor Td, Q4.1D and ( |4.2| ) can be rewritten as 

y=p(t)-y[l + rf l (y-z) z ] 
e z = -{z - 1) - e z [l + rf{y - z) 2 ] , 

where £o = T r/ T d ; V 2 = r d(«T#) 2 9/U, and p(t) is proportional to the freshwater flux 
F(t), with a factor olsSqt&/ (oltQH). Cessi uses the estimates rj 2 ~ 7.5, r r ~ 25 days and 
Td ~ 219 years. This yields Eq ~ 3 x 10 -4 , implying that (fOj) is a slow-fast system. 



Tihonov's theorem [36] allows us to reduce the dynamics to the attracting slow manifold 



z = 1 + 0(eq). To leading order, we thus find 

y = -y[l + r ] 2 (y-l) 2 ] +p(t). (4.6) 

Stochasticity shows up in this model through the weather-dependent term p{t). To 
model long-scale variations in the typical weather, we will assume that p(t) can be repre- 
sented as the sum of a periodic term p(t) and white noise, where the period 1/e of p(t) is 
much longer than the diffusion time, which equals 1. (Recall that we have already rescaled 
time by a factor of r^.) We thus obtain the SDE 

dy t = f(y t ,t)dt + a dW t , where f(y,t) = -y[l + r, 2 {y - l) 2 ] +p{t). (4.7) 

Note that / has an inflection point at y = 2/3, and that 



P(*) "3-^7^ 



+ 



3 ' 



(4.8) 



which derives from the Ginzburg-Landau potential ( |1.8| ) with parameters /u = {rf /3 — 1) 
and A(t) = T](p(t) — 2/3 — 2rj 2 /27). As we already know, the potential has two wells if and 
only if A 2 < A 2 = 4/i 3 /27, which means, for rj 2 = 7.5, that p G [0.96, 1.48]. The double-well 
potential is symmetric for p = po = 2/3 + 2r/ 2 /27 ~ 1.22. 

For a deterministic forcing given by X(t) = K cos2iret, the SDE for x = r](y — 2/3) 
becomes, on the time scale 1/e, 

dx t = -[px - x 3 + K cos 2vrt] di + dW t , (4.9) 
e ye 

where cr = do?/. This SDE is of the same form as (|3.8j). While in Section |3[ we assumed 
ii' < A c , we will now allow X to exceed A c , so that the difference ao = K — A c may change 
sign. (Note that in Section 0, ao had the opposite sign.) 



In the deterministic case a = 0, Equation (|4.9| ) has been used to model a laser [23|, and 
a similar equation describes the dynamics of a mean-field Curie- Weiss ferromagnet [37]. 
In the limit of infinitely slow forcing, solutions always remain in the same potential well 
if K < A c . If K > A c , however, the well tracked by x% disappears in a saddle-node 
bifurcation when \\{t)\ crosses A c from below, causing xt to jump to the other well, which 
leads to hysteresis, see Figure 

For positive e, the system does not react immediately to changes in the potential, so 
that the hysteresis cycles are deformed. One can show p3|, ITS that 
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K u w > J 



FIGURE 4. The potential V(a;,i) = jar - |ar - X(t)x, with A(i) = -ftTcos(27rt), when K 
exceeds the threshold A c . In the deterministic case, with s <C 1, the overdamped particle 
jumps to a new well whenever |A(t)| becomes larger than A c , leading to hysteresis. Larger 
values of s increase the size of hysteresis cycles, but additive noise of sufficient intensity 
decreases the size of typical cycles, because it advances transitions to the deeper well. 

For K ^ A c + 0(e), xt always tracks the same potential well, at a distance at most of 
order e/ y^Jao] if oo ^ — e, and of order i/e if |oo| is of order e. 

For K ^ A c + 0(e), xt is attracted by a hysteresis cycle, which is larger than the static 
hysteresis cycle; in particular, xt crosses the A-axis when \(t) = K cos 2nt = A , where 
A satisfies 



|A l 



with ao = K — A c 



(4.10) 



Additive noise will also influence the shape of hysteresis cycles, because it can kick the 
state over the potential barrier, as has been noted in |28| in the context of the thermohaline 
circulation. For positive noise intensities a, the value A at which x% crosses the A-axis, 
becomes a random variable. Assume for instance that we start at time to = 1/4 in the 
right-hand potential well. We define 



H = inf {t G [± U : x t (u) < o}, A°(w) = X(r°(u)), 



(4.11) 



with the convention that t°(uj) = oo and X°(u>) = oo if Xt(uj) > for all t G [4, f]- We 
thus have r° G [|, |] U {00} and A G [— K,K] U {00}. We will indicate the parameter- 
dependence by A = A°(e, cr), keeping in mind that this random variable also depends 
on ao and /i. In the deterministic case, A°(e,0) = 00 if K ^ A c + 0(e), and A°(e, 0) 
satisfies flOoD if if ^ A c + 0(e). 

As we know from the previous section, for K < A c , there is an amplitude-dependent 
threshold noise level a c such that during one period, x% is unlikely to cross the potential 
barrier for a <C cr c , while it is likely to cross it for a 3> cr c . In fact, in the latter case, there 
is a large probability to cross the barrier a time of order <r 2 / 3 before the instant t = 1/2 
of minimal barrier height, when A is of order A c — <t 4//3 . In that case, the hysteresis cycle 
will be smaller than the static cycle. A similar distinction between a small-noise and a 
large-noise regime exists for large-amplitude forcing. 

It turns out that the distribution of A can be of three different types, depending on 
the values of the parameters (c. f. Figure [| and Figure |6|): 
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Figure 5. Typical random hysteresis "cycles" in the three parameter regimes, (a) Case I: 
Driving amplitude K and noise intensity a are too small to allow the path to switch poten- 
tial wells, (b) Case II: For large amplitude but weak noise, the path tracks the deterministic 
hysteresis cycle, which is larger than the static one. (c) Case III: For sufficiently strong 
noise, the path can overcome the potential barrier, so that typical hysteresis cycles are 
smaller than the static one. 



• Case I — Small-amplitude regime: ao ^ const e and a ^ (|oo| Ve) 3 / 4 . 

Then xt is unlikely to cross the potential barrier, and there are constants C, k > 
such that (see []|, Theorem 2.3]) 

P{A° < OO} < - e -<\a \Vef/ya* _ (412) 

The probability to observe a "macroscopic" hysteresis cycle is very small, as most 
paths are concentrated in a small neighbourhood of the bottom of the right-hand 
potential well (Figure ||a). 

• Case II — Large-amplitude regime: ao ^ const e and a ^ (e^/oq) 1 ^ 2 . 

This regime is actually the most difficult to study, since the deterministic solution 
jumps when \X(t)\ — A c X (e^fa® ) 2 / 3 , and crosses a zone of instability before reaching 
the left-hand potential well. One can show, however, that |A°| is concentrated in an 
interval of length of order (e^fa®) 2 ^ around the deterministic value |9|, Theorem 2.4]. 
More precisely, there are constants C, k > such that 

P{|A°| < A c -L} < £ e -^ 3/2 /- 2 (4.13) 

for (e^) 2/3 < L < L /|logK/^)l, and 

F{|A°| < Ac + L^e^) 2 / 3 } s$ £ e - K£ v^/- 2 5 (4.14) 

where the constants Lq,Li > are independent of the small parameters. Hence it 
is unlikely to observe a substantially smaller value of |A°| than the deterministic one, 
provided a <C {e^/oq) 1 ^ 2 . On the other hand, there is a constant L2 > L\ such that 

P{|A°| > A c + L) 3e -^/(^(^) 2/3 |i°g(^)l) (4.15) 

for all L > L 2 (e v / ^ u ') 2/3 . As a consequence, the vast majority of hysteresis cycles will 
look very similar to the deterministic ones, which are slightly larger than the static 
hysteresis cycle (Figure ||b). 
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FIGURE 6. The three hysteresis regimes, shown in the plane driving-amplitude-noise- 
intensity, for fixed driving frequency. The insets sketch the distribution of the random value 
A of the forcing A(i) when Xt changes sign for the first time. In Case I, such transitions are 
unlikely. In Case II, |A°| is concentrated in an interval [A c +ii(e v / ao ) 2 / 3 , Ac + T^e-^/ao) 2 / 3 ] 
containing the deterministic value |A°(e,0)|. The broken curve indicates that we do not 
control the distribution inside this interval. In Case III, |A°| is concentrated around a value 
A which is smaller than A c by an amount of order <r 4 / 3 . The distribution decays faster to 
the right, with a width of order e (actually, e|logtr|) if ao ^ e or a ^ |ao| 3 ^ 4 , and of order 
e^/a^/a 2 / 3 if a ^ e and a < |a | 3 / 4 . 



Case III — Large- noise regime: Either ao ^ e and a ^ (|ao| V tr) 3 / 4 or ao ^ £ and 

In this case, the noise is sufficiently strong to drive xt over the potential barrier, 
with large probability, some time before the barrier is lowest or vanishes, leading to a 
smaller hysteresis cycle than in the deterministic case (Figure ||c). It turns out that 
I A 1 is always concentrated around a (deterministic) value A satisfying A c -Ax cr 4 / 3 . 
It follows from |9|, Proposition 5.1] that 

P{|A°| < A - L) < ^ 6 - kL3/2 / ct2 +-e— 4/3 /(e|iog-|) (416) 
s 2 



for ^ L C A and 



F{|A°| > A + L) ^ - e -^/(e\loga\) ( 41 ^ 



for positive L up to K — A if ao ^ e. If ao e, the same bound holds for L ^ 
A c — A, while the behaviour for larger L is described by ( 4.15| ). The estimates ( |4. 16| ) 
and ( 4.1?1 ) hold if ao ^ e or a > a^ 4 . In the other case, two exponents are modi- 
fied: <7 4//3 /(cT|log o" |) is replaced by a 2 /{e^/liQ\\oga\), and L/(e|log <r|) is replaced by 
a 2 / 3 L/( e ya7J|logcT|). 

Note that in all cases, the distribution of A decays faster to the right than to the 
left of A, and it is unlikely to observe A larger than A c , except when approaching the 
lower boundary of Region III. 
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In some physical applications, for instance in ferromagnets, the area enclosed by hys- 
teresis cycles represents the energy dissipation per period. The distribution of the random 
hysteresis area can also be described, and bounds on its expectation and variance can be 
obtained. We refer to || and []|, Section 4] for details. 

For Stommel's box model, the above properties have two important consequences. 
First, noise can drive the system from one stable equilibrium to the other before the 
potential barrier between them disappears, so that a smaller deviation from the mean 
freshwater flux than expected from the deterministic analysis can switch the system's 
state. Second, this early switching to the other state is likely only if the noise intensity 
exceeds a threshold value (which is lowest when the amplitude K is close to A c ). Still, 
the system spends roughly half of the time per period in metastable equilibrium in the 
shallower well. 



5 Delay 

Convective motions in the atmosphere can be simulated in a laboratory experiment known 
as Rayleigh-Benard convection. A fluid contained between two horizontal plates is heated 
from below. For low heating, the fluid remains at rest. Above a threshold, stationary 
convection rolls develop. With increasing energy supply, the angular velocity of the rolls 
becomes time-dependent, first periodically, and then, after a sequence of bifurcations de- 
pending on the geometry of the set-up, chaotic. For still stronger heating, the convection 
rolls are destroyed and the dynamics becomes turbulent. 

Lorenz' famous model |2(| uses a three- modes Galerkin approximation of the hydro- 
dynamic equations. The amplitudes of these modes obey the ODEs 

X = Pr(y-X) 

Y = rX -Y - XZ (5.1) 
Z = -bZ + XY. 

Here X measures the angular velocity of convection rolls, while Y and Z parametrize the 
temperature field. The Prandtl number Pr > is a characteristic of the fluid, b depends 
on the geometry of the container, and r is proportional to the heating. 

For ^ r ^ 1, the origin (X,Y,Z) = (0,0,0) is a global attractor of the system, 
corresponding to the fluid at rest. At r = 1, this state becomes unstable in a pitchfork 
bifurcation. Two new stable equilibrium branches (±y / 6(r — 1), ±-v/6(r — 1), r — 1) are 
created, which correspond to convection rolls with the two possible directions of rotation. 
We will focus on this simplest bifurcation, ignoring all the other sequences of bifurcations 
ultimately leading to a strange attractor (see for instance f32fl ). 

We are interested in the situation where r = r(et) grows monotonously through r(0) = 
1 with low speed e (e. g. r = 1+et). Near the bifurcation point, one can reduce the system 
to an invariant center manifold, on which the dynamics is governed (c. f. Jl0|] ), after scaling 
time by a factor e, by the one-dimensional equation 

e^L = fi(t)x + c(t)x 3 + 0(x 5 ). (5.2) 
at 

Here p(t) = a(t) + 0{e), where a{t) = § [-(Pr + 1) + ^(Pr + l) 2 + 4Pr(r(t) - 1)] is the 
largest eigenvalue of the linearization of (|5.1|) at 0, which has the same sign as r(t) — 1, 



and c(t) is negative and bounded away from zero. The right-hand side of (5.2) derives 



18 



Figure 7. The potential y(x,t) = jx 4 — ^u(t)x 2 transforms, as \x changes from negative 
to positive, from a single-well to a double-well potential. In the deterministic case, an 
overdamped particle stays close to the saddle for a macroscopic time before falling into 
one of the wells. Noise tends to reduce this delay. 



from a potential similar to the Ginzburg-Landau potential (1.8) with A = 0, which remains 
symmetric while transforming from a single-well to a double-well as fj,(t) becomes positive, 
see Figure |7[ 

The solution of fl5,2| ) with initial condition xq > for to < can be written in the 
form 



x t = ¥>(zo,t)e Q(t, '° )/e , «(Mo)= / ds, (5.3) 



to 

with < (f(xQ,t) sj xo for all t. Thus xt is exponentially small if a(t,to) is negative. The 
important point to note is that a(t,to) can be negative even when a(t) is positive. For 
instance, if fx(s) = s, then a(t, to) = \{t 2 — t^) is negative for to < t < —to- Thus x t 
will remain exponentially close to the saddle at x = up to time —to after crossing the 
bifurcation point. This phenomenon is called bifurcation delay. It means that when r is 
slowly increased, convection rolls will not appear at r = 1, as expected from the static 
analysis, but only for some larger value of r, which depends on the initial condition. 

It is clear that the existence of a delay depends crucially on the fact that Xt can 
approach the saddle exponentially closely, where the repulsion is very small. Noise present 
in the system will help kicking xt away from the saddle, and thus reduce the delay. The 
question is to determine how the delay depends on the noise intensity a. 

For brevity, we will illustrate the results in the particular case of a Ginzburg-Landau 
potential, with dynamics governed by the SDE 

dx t = -[fx(t)x t -xf]dt+-^=dW t . (5.4) 
e ye 



The case without the term —x\ has been analysed by several authors [35], 22], with 

the result that the typical bifurcation delay in the presence of noise behaves like yjlogoj. 
The results in || cover more general nonlinearities than — x 3 . 

We assume that \x{€) is increasing, and satisfies /z(0) = 0, //(0) const > 0. For 
simplicity, we consider first the case where xj starts at a time to < at the origin x = 0. 
From the results of Section we expect the paths to remain concentrated, for some time, 



in a set whose width is related to the linearization of (5.4) around x = 0. We define the 
function 

v(t) = v e 2a{t),£ +— f e 2a{t ' s)/£ ds, where a(t) = a(t, 0). (5.5) 

£ Jtn 
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For a suitably chosen vq x a 2 /\fi(to)\, one can show that v(t) is increasing and satisfies 
v(t) 



'a 2 /\n{t)\ for t < t < -y/i 

a 2 /y/e~ for -y/i ^t^y/i (5.6) 

a 2 e 2*(t)/e fort^Ve- 



Note that although the curvature of the potential at the origin vanishes at time 0, 

v(t) grows slowly until time y/e after the bifurcation point, and only then it starts growing 
faster and faster. 

We now introduce, as in Section 0, the set 



B{h) = {(x,t): \x\ < hy/€(jj}. (5.7) 

Then one can show (see H Theorem 2.10]) the existence of a constant ho > such that 
the first-exit time tb(^) of xj from B(h) satisfies 

f{r B{h) <t}^C(t,e)e- Kh2 (5.8) 

for all h sj hoa/v(t), where 

C(t,e) = ^yW)|ds + oQ, and k = ^ - 0^) - (^^-j- (5-9) 

The paths are concentrated in B(h), provided h^a 2 /v(t) 2 ^ h 2 3> log C(t,e). As a conse- 
quence, we can distinguish between three regimes, depending on noise intensity: 

• Regime I: a ^ e~ K / £ for some K > 0. 

The paths are concentrated near x = at least as long as 2a(t) <C K. This implies 
that there is still a macroscopic bifurcation delay. 

• Regime II: e _1 / £P ^ a <C fo r some p < 1. 

The paths are concentrated near x = at least up to time y/e, with a typical spreading 
growing like a/\/\fi(t)\ for i —y/z-, and remaining of order a/e 1 ^ for |t| ^/e. 

• Regime III: a ^ ^/i. 

The paths are concentrated near x = at least up to time — a, with a typical spreading 
growing like aj ■v/|/x(i)|. Near i = 0, the potential becomes too flat to counteract the 
diffusion, and as t grows further, paths keep switching back and forth between the 
wells, before ultimately settling for a well. 

Similar results hold if xt starts, at to < 0, away from x = 0, say in xo > 0. Then the set 
B(h) is centred at the deterministic solution x^ et (with the same initial condition), which 
jumps to the right-hand well when a(t,to) becomes positive, see Figure ||. In Regime I, 
with K sufficiently large, the majority of paths follow xf et into the right-hand potential 
well. 

It remains to understand the behaviour after time y/e in Regime II. To this end, we 
introduce the set 



V(g) = {(x,t):t^^,\x\ < y/(l-g)n(t)}, (5.10) 

depending on a parameter g £ [0,2/3). The set T>(0) contains the points lying between 
the two stable equilibrium branches ±yj /x(t). One can show (see [||, Theorem 2.11]) that 
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Figure 8. A sample path x t of the SDE ( pT4| ) with fi(t) = t, for e = 0.01 and a = 0.015 



The deterministic solution xf , starting in xq > at time to, jumps to the right-hand 
well, located at x*(t) = Vt, at time |io|. Typical paths stay in the set B(h), whose width 
increases like ha/{yJ\t\\J e 1 ^ 4 ), until time ^fe after the bifurcation. They leave the domain 



T>(g) (shown for g = 2/3) at a random time r = Tj}/ g \, which is typically of order ^/e|logcr|. 
After leaving T>(g), each path is likely to stay in a strip of width of order ha/y/t, centred 
at a deterministic solution approaching either +x*(t) or —x*(t). 

if g € (0,2/3) and <7|logcr| 3 / 2 = 0(y/e), then the first-exit time Tx>r g ) of xt from T>(g) 
satisfies 

Hr v{Q) >t}^ C(t, e) 1 8 1 , (5.11) 

CT \J\ — e -2ga(t,^/e)/e 

where 

C(t,e) = const M (t)(l + ^i^£T). (5.12) 



The estimate ( |5.11| ) shows that paths are unlikely to stay in T>(g) as soon as t satisfies 
ga(t,y/e) 3> e|logcr|. Since a is quadratic in t, most paths will have left T>(g) for 



t » Vepogoi. (5.13) 

Once xt has left T>(g), one can further show that it is likely to track a deterministic 
solution which approaches the bottom of one of the potential wells. Assume for instance 
that xt leaves T>(g) through the upper boundary, at a random time r = Tx>( e )- Then, for 
1/2 < g < 2/3, |], Theorem 2.12] shows that the deterministic solution xf et,T , starting in 
the same point at time r, approaches the bottom of the well at \J fx(t) like e//i(t) 3//2 + 
\J (j>(t) e~ r?a (*' T )/ £ , where r] = 2 — 3g, and the path xt is likely to stay in a strip of width 
a I \J fi(t) around x dct,T . Thus after another time span of the form ( p. 13 ), most paths will 
have concentrated near the bottom of a potential well again. 

We note that different kinds of metastability play a role here. First, paths remain 
concentrated for some time near the unstable saddle. Second, they will concentrate again 
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near one of the potential wells after some time. Some paths will choose the left-hand well 
and others the right-hand well (with probability exponentially close to 1/2 in Regime II), 
but all the paths which choose a given potential well are unlikely to cross the barrier again. 
In fact, one can show that if fi(t) grows at least linearly, then the probability ever to cross 
the saddle again is of order e~ const l a . If we start the system at a positive to in one of the 
wells, the distribution will never approach a symmetric bimodal one. 

In the case of the Rayleigh-Benard convection with slowly growing heat supply r(et) 
and additive noise, these results mean that exponentially weak noise will not prevent the 
delayed appearance of convection rolls. For moderate noise intensity, rolls will appear 
after a delay of order y/\\oga\/e, which is considerably shorter than the delay in the 
deterministic case which is of order 1/e. The direction of rotation is unlikely to change 
after another time span of that order. For strong noise, convection rolls may appear early, 
but their angular velocity will fluctuate around zero until a time of order a/e after the 
bifurcation before settling for a sign, and even then occasional changes of rotation direction 
are possible. 
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